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The decay of a metastable system is described by extending Kramers' method to the 
quantal regime. For temperatures above twice the crossover value we recover the result 
known from applying Euclidean path integrals to solvable models. Our derivation is not 
restricted to a linearly coupled heat bath of oscillators, and thus applicable to nuclear 
systems. 
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In his famous paper Kramers [1] derived a formula which took into account the im- 
plications of dynamics on the decay rate of a metastable system. Applied to the situation 
portrayed in Fig.l it can be written as: 

^ = ^( V / TT^-^e-^ (1) 

Here, (3 = (ZcT) -1 , B is the barrier height and w a the frequency of the motion around the 
potential minimum (fulfilling the relation M a w 2 = C a = d 2 V(Q)/dQ 2 \ a with M a being 
the inertia). The rj b = i h /{2y/M b \ C b |) measures the effective damping rate, calculated 
at the barrier from the friction coefficient 75, the inertia M b and the stiffness coefficient 
C b = d 2 V(Q)/dQ 2 \ b . 

Expression (1) is valid for damping rates r\ b not smaller than about 0.2. Furthermore, 
one has to assume the barrier big enough to keep the decay rate R sufficiently small. Under 
such conditions the process can be viewed as quasi-stationary with the outward flux j b at 
Q b being constant for macroscopically large times. The decay rate R is then simply given 
by the ratio of the current divided by the number of particles N a caught in the well near 

Qa, R = jb/N a . 

It has always been a challenge to extend Kramers' result to the quantal regime. How- 
ever, it was only after one had learned to apply the instanton trick to this problem that a 
decent solution could be found. Among the vast literature on this subject of "dissipative 
tunneling" we would only like to refer to the classic paper [2], to [3] as a general recent 
review and to [4] as a reference to formal details. Common to all these approaches is 
the application of the technique of path integrals for imaginary time propagation. This is 
feasible if for the basic Hamiltonian one assumes a form like: 

H( Xi , Pi , Q, P) = H B (x llPt ) + H SB {xi,Pi, Q, P) + H S (Q, P) (2) 

in which two restrictions are imposed. First, the coupling is assumed to be of the form 
Hsb = QF with F = J2i c i x i- Secondly, the heat bath must be represented by a set of 
oscillators. It is needless to stress that some measures have to be taken in order to render 
the final equations of motion irreversible. 

A first major step to unify these results with the usual concepts of transport theory 
has been undertaken in [5, 6]. Still using path integrals, the author has been able to 
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obtain from real time propagation a quantal version of Kramers' stationary solution for 
the inverted oscillator, which he could then use to rederive the form of the known quantal 
correction factor. This derivation is valid above a critical temperature T c = 2T with T 
being the so called crossover value [3, 4]. We will come back to the physical significance 
of this restriction. 

A full-fledged transport theory has been formulated in [7] within a quasi-classical 
approach based on a quantal Langevin equation. Unfortunately, the authors have not 
been able to establish a direct connection to the results mentioned above. This will be 
the main goal of the discussion to come. We will proceed by exhibiting first the formal 
details, and then discuss later both the physical conditions of the derivation as well as 
consequences for possible applications. 

Let us begin by formulating the decay rate in terms of Wigner functions k w ((3; Q, P). 
Borrowing from the ideas of Kramers we may write: 



Here k^ oh (Q, P) is meant to represent the full global solution to our problem. Under the 
conditions mentioned already above, and for not too small temperatures (see below), this 
global function can be approximated by local ones valid at the barrier or at the potential 
minimum, respectively. They are obtained by appropriate linearizations. Trivially an 
overall normalization factor drops out of eq.(3) and it is almost obvious that the local 
solutions should not be normalized in themselves. In this context it is useful, for a moment, 
to think of a stationary situation. Then kg[ oh (Q, P) would simply be the Wigner transform 
of fcgiob = exp(— (3H). As a consequence of this transformation one gets a prefactor (to the 
exponential factor depending on the coordinates and momenta), which in the end may be 
found from the partition function 2 g i b- The linearization procedure we spoke of will in 
such a case have to be performed on Ti. In this way it will lead to both a local partition 
function as well as to a local prefactor Q. Both will depend on the local frequency rendering 
obvious that the Q will have to be different at the well and at the barrier. Our situation 
differs from the one just described only in the sense that we have to deal with a dynamical 
situation, albeit in quasi stationary approximation, whose solutions must be found from 
transport equations. 




(3) 
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In [8] a Hamiltonian of type (2) has been used to derive a transport equation for 
local collective motion. This Hamiltonian can be obtained self-consistently by applying 
an appropriate linearization procedure with respect to the macroscopic degrees of freedom 
(see [9]). No restrictions need to be imposed on the "bath" Hamiltonian and the F 
appearing as a factor in Hqb (see text below eq.(2)). As we shall see, for the present 
purpose it suffices to expand around values Qo which correspond to extremal points of the 
effective potential (for the general case see eg. [8] or [10]). Then the result can be written 
as: 



|/(,,p,«) 



d P d „ d P „ d 2 „ d 
+ 1HS C <1 + + D qp —- + D : 



dqM dP dPM ' qp dqdP pp OPdP 



f{q,P,t) (4) 



with q = Q — Qo. (With respect to P we assume from the start that our Hamiltonian 
is at most quadratic). Besides the transport coefficients for average (local) motion, there 
appear the two diffusion coefficients D pp = {l/M)Y, pp and D qp = CJ2 qq - 0-/ M )Y, pp i 
with the equilibrium fluctuation J2 qq given by the fluctuation dissipation theorem (FDT): 

= St coth (£) *» < 5 » 

Similar expressions hold for J2 PP an d J2 qp ^ f° r each appearance of the kinetic momentum 
P a factor Mu is multiplied to the integrand. As a consequence, the J2 qp vanishes because 
of a simple symmetry, which is the reason for D qq = 0. 

The X qq (uj) is the dissipative part of the susceptibility Xqqi 00 ) which measures the 
response of (q) u to an external "field" q e xt(<^) i n linear order. These functions contain 
information about all possible modes of our system, as given by the Hamiltonian (2). 
Clearly, we are interested only in the collective ones, namely those whose average behaviour 
we have parametrized in terms of the transport coefficients M, 7 and C, defined for local 
motion. To render the dynamics of the fluctuations consistent, we have to calculate the 
diffusion coefficients from the corresponding approximations. That means to replace in (5) 
the locally valid x qq (uj) by the dissipative part of the oscillator response function being 
defined by: Xosc( UJ ) c l( UJ ) = ( — ^ 2 M — iwy + C)q(u) = —q ext (uj). We are aware that this 
restriction will finally cause the expected problem of a diverging Y1 PP - To regularize the 
integral appearing there one has to introduce one other constant. For example, this could 



be done by way of a cut off [8] or through a Drude regularization (cf.fll] and [12]). For 
our present purpose this problem is not really relevant. 

We are ready now to evaluate the decay rate. Let us look first at the denominator of 
(3). The linearization we spoke of before means to put kg[ oh (Q ~ Q a , P) ~ k^ c (q, P) and 
to use the following stationary solution of (4): 

Cc(^) = Ca^ fl exp(--^- - ^=^) (6) 



with the prefactor given as 2n( a = Z a / \jYT qq YTppi which normalizes (6) to Z a as dis- 
cussed above. (Here and in the following the Z : s always refer to the effective macroscopic 
part of the partition function). In order to calculate the number of particles N a we may 



replace A by oo, provided we have \Q a \ ^ yYlqq- Next we turn to the barrier region. 
Here the stiffness is negative and the basic expressions for the diffusion coefficients be- 
come somewhat delicate, as the "equilibrium fluctuations" loose their immediate physical 
meaning. Nevertheless, all formulas can be continued analytically by replacing 



everywhere. In this way the J2 qq becomes negative, as already known from Kramers' case 
of high temperatures. For the prefactor at Qb this implies to write 2n(b = Zbj yYlqqYlpp = 



As for the functional dependence on Q and P, a form like the one obtained from (6) 
still solves our equation of motion (4), now written for the local motion around Qb = 0. 
However, this would not be a decent approximation to k^( oh (Q,P). Indeed, since the 
barrier is unbound to the right we expect a solution with a finite current. Such a solution 
has been displayed in [13] and applied to computations of interest in nuclear fission, but not 
to the decay rate. Its form had been found before in [5] within a path integral description. 
In our present notation it reads 

q 2 \ f P ~ Aq du - - 2 



k$ oh (Q ~ Qb,P) - k?(q,P) = C & exp(-— - - « ) / " exp[-|-] (8) 



jpp ^,qq 



with a = -Y.'qqA 2 - El and M b w b (y/T+vl - = Ep P /("E^) 



The final expression for the decay rate is easily obtained. After putting into (3) both 
(6) and (8) the integrals can be calculated by elementary methods to give: 



(9) 



As observed already in [5], the quantum correction factor {q can be expressed by the 
absolute value of the ratio of the local partition functions at a and b. In the remaining 
part of our letter we will discuss the physical relevance of our result together with a critique 
of the present derivation. 

1) Comparison with previous derivations. The partition functions appearing in (9) are 
readily calculated for the model of the bilinear coupling to a set of oscillators, for example 
with path integrals (cf. e.g. [4] and our discussion below eq.(2)). Again take first the case 
of the potential minimum where one gets: 

-1/2 

Za=M\ I! (10) 




with the eigenvalues A^ = M a v\ + \v n \la(Wn\) + M a w^. Here, v n = 2nn/hf3 and Af is 
the usual reference factor being independent of w a ,M a and 7 a . (We follow the notation 
of [4] allowing frequency dependent friction coefficients 7(0;) for an eventual Drude regu- 
larization). The Z\> is obtained from the transformation (7). We observe that the mode 
n = has a negative eigenvalue, which for this model proves the conjecture about Z^ 
stated above, namely that Z^ becomes purely imaginary. For temperatures above Tq, with 
the corresponding /? being a solution of A^ = X b _ 1 = 0, all other eigenvalues are positive. 
For the quantum correction factor one thus obtains: 

nii^ + ^K)-^ 2 [n) 

with T = 7/M. This expression has been derived before within the formulation of "dis- 
sipative tunneling", but for the case M a — Mb and 7 a = 7*, *. We present it in Fig. 2 
for a whole range of effective damping factors rjb- At the crossover temperature Tq this 



* The case of variable inertia and friction has been treated in [14] within a phenomeno- 
logical ansatz; see our discussion to come in 4). 
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formula diverges, but it can be regularized after considering deviations of third order from 
the saddle points [4]. As we shall see now our derivation is restricted to temperatures of 
approximately twice that value, T > 2T . 

2) The limits of our approach. In our discussion above we have been very sloppy 
about any convergence problem appearing in the integrals. Looking back at formula (3) we 
realize that problems may come from the numerator of this expression. We certainly need a 
positive Ylpp- However, the more crucial quantity turns out to be the Y^ qq - ^ n Kramers'case 
it is known to be negative such that the Gaussian factor in (8) increases with increasing 
q. To study these expressions in general and for sizable friction is somewhat elaborate, 
at least for the regime where quantal effects are important. For a first orientation, let us 
therefore look at weakly damped motion, to postpone a discussion of the general case to 
future publications (see e.g. [10]). Evaluating the Ylnu f rom (5) to zero order in 7 one gets: 



which is nothing else but the quantal virial theorem continued analytically to the imaginary 
frequency of the barrier, making apparent the fact that the "equilibrium fluctuation" in 
q is negative at the barrier. (To get Kramers' case of high temperature we just have to 
replace T* by T). For a situation like the one behind (12) the condition we have to pose 
clearly is T > T c = 2T = hwb/n. It guarantees that J2 PP * s positive and J2 qq negative. 
Notice, that for this situation it is easy to observe that also the a of eq.(8) is positive. At 
T = T c = 2T Q the Gaussian factor appearing in (8) becomes indifferent as function of q 
and regularization procedures would have to be employed [15,6]. 

3) On the quantum nature of this approach. Our result clearly demonstrates that 
the transport equation (4) contains quantum effects. It is clear that the latter are there 
on a semiclassical level only, as global motion is treated in a locally harmonic fashion. 
Nevertheless it is interesting to see the effects contained in a differential transport equation. 

Without any doubt this equation causes problems for the basic uncertainty relations 
at short times. This is no surprise as we have left out memory effects as well as the 
implications of initial correlations. That both these effects must be there can be studied 
in the oscillator model again. One may look, for instance, at the results of [16] for the 
time evolution of the second moments, say for the case of a positive stiffness. Compared 




(12) 
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with the evolution obtained from (4), one sees that both agree for large times. This is 
to be expected, of course, from the very construction of the transport equation [8]: The 
diffusion coefficients are chosen such as to guarantee the correct equilibrium, as defined by 
the FDT. 

The lesson we learn from our present results is that these features prevail also for 
the case of the inverted oscillator. Again, this does not come unexpected. One knows 
both from [17] and [5] that for the inverted oscillator time dependent solutions "relax" 
to stationary ones. Apparently, in the regime of temperatures we look upon here, it is 
the long time behaviour of the solution which matters — and this behaviour contains the 
adequate quantum feature. Further details will be discussed elsewhere, see e.g. [10]. 

4) Advantages of our approach. Although our method does not apply at very low 
temperatures, we hope to contribute to the understanding of quantum effects in dissipative 
dynamics. We see at least the following three favorable circumstances. 

i) It should be of some theoretical interest to see effects of "dissipative tunneling" come 
out of a transport equation proper. 

ii) Our derivation is not based on the assumption of a linear coupling to a heat bath which 
does not change with Q. It is commonly understood that such a model would be quite 
unrealistic for nuclear fission, a case discussed already by Kramers. In a Hamiltonian 
like (2) the Xi,pi represent the dynamics of the nucleons. For a first approximation there 
is no way around a mean field approach, which poses important constraints both on the 
Hamiltonian itself as well as on the treatment of the coupling (see [9] and further literature 
cited there). First of all, Hsb must have a term involving the collective momentum P. 
Secondly, the dependence of Hb and Hqb on the Xi cannot just be of first and second order. 
This feature becomes already evident by looking at the simple but well examined case of 
multipole vibrations around a potential minimum. It is only for dipole excitations that the 
Hqb can be taken linear in the X{. Thirdly, for large scale motion both Hb as well as Hsb 
must vary along the collective path. Within the locally harmonic approach one is able to 
incorporate these effects both for the construction of the Hamiltonian 7i(xi,Pi, Q, P) itself 
[9] as well as for the derivation of the equations of motion [8]. 

We have mentioned already ref.[14]. There the effects of variable inertia and friction on 
the decay rate have been taken into account by the following procedure: It was assumed 
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that a Hamiltonian like (2) exists for global motion, with a factorized form for the Hsb', the 
latter was supposed just to depend on the coordinates, thus neglecting self-consistency; 
some phenomenological ansatz was made both for the form factor of Hsb as well as for 
the (unperturbed) inertia in Hs, neglecting the renormalization of the conservative forces 
through the coupling; the collective variables were linearly coupled to a heat bath of 
oscillators with fixed frequencies. 

iii) Our approach may allow for interesting practical applications, even beyond the mere 
calculation of decay rates. For instance, one is not bound to look at cases where the 
temperature stays constant. Moreover, as we look at real time propagation our transport 
equation can easily be used for Monte Carlo simulations. Notice, that our locally harmonic 
approximation concurs ideally with the small time steps employed in these calculations. 
Such simulations have proven successfully in many fields of physics. Let us only allude to 
the following examples where instabilities play a dominant role: once more nuclear fission, 
chemical and surface reactions, nucleation and spinodal decompositions. The lesson we 
learn from the present study is how simple it can be to incorporate quantum effects: 
Suppose one wants to treat time evolution of the average values on the basis of classical 
physics. The use of proper diffusion coefficients restores quantum features on the level 
of the second moments, a procedure which is quite easy to perform within the locally 
harmonic approximation. In the case discussed here average dynamics is represented by 
our collective variable Q(t), but an extension to time evolution of the general mean field 
is straightforward. For the latter the appropriate classical means is the Landau- Vlasov 
equation. A linearization allows to deduce response functions which in turn determine 
quantum statistical fluctuations and hence generalized diffusion coefficients (cf.[18]). 

Acknowledgments: We would like to thank H. Grabert, P. Frobrich and R. Sollacher for 
fruitful discussions and valuable suggestions. 
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FIGURE CAPTIONS 



Fig.l: A schematic plot of the effective potential. 

Fig. 2: The quantum correction factor {q for the case of both constant inertia M a = Mb and 
friction 7 a = 75, and the following values of the effective damping rate r\b = 10, 5, 2, 1, 0.5, 
from left to right, respectively. The arrows show the corresponding crossover temperatures 
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